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SECTION  I 


INTRODUCTION 


Bounded-wave  electromagnetic-pulse  (EMP)  simulators  provide  one  means  to 
generate  simulated  EMP's  for  testing  the  EMP  responses  of  aircraft  and  missiles 
on  the  system  level  (ref.  1).  A  simulator  can  be  designed  to  launch  an  electro¬ 
magnetic  pulse  into  the  simulator  volume  which  closely  approximates  the  actual 
nuclear  EMP.  This  ability,  however,  does  not  imply  that  the  actual  interaction 
between  the  nuclear  EMP  and  the  test  object  is  also  thereby  automatically  closely 
simulated.  In  the  simulator  there  exists  an  additional  interaction  between  the 
simulator  structure  and  the  test  object  which  is  absent  in  the  actual  EMP 
encounter  in  free  space.  This  interaction  adversely  affects  the  quality  of 
the  simulation.  In  order  to  interpret  and  utilize  correctly  the  simulator 
test  data,  it  is  necessary  to  estimate  and  make  allowance  for  this  simulator/ 
test-object  interaction. 

This  report  presents  an  analytical  study  of  the  simulator /test-object 
interaction.  It  employs  an  idealized  model  with  simplified  geometries  to  keep 
the  mathematics  manageable.  The  bounded-wave  EMP  simulator  is  modeled  by  a 
pair  of  parallel,  perfectly  conducting,  flat  plates  of  infinite  extent,  as 
depicted  in  figure  1.  Such  a  structure  is  an  electromagnetic  waveguide.  Let 
the  spacing  between  the  two  plates  be  denoted  by  h  .  A  rectangular  coordinate 
system  can  be  set  up  such  that  the  lower  plate  is  at  z  ■  0  and  the  upper  plate 
at  z»h  .  The  test  object  is  taken  to  be  a  hollow  cylinder  or  tube  of  cross- 
sectional  radius  a  and  length  L  .  It  is  considered  to  be  perfectly  conducting 
and  of  infinitesimal  thickness.  The  axis  of  the  tube  is  oriented  perpendicular 
to  the  two  plates,  and  can  be  taken  to  coincide  with  the  z-axis.  The  length 
L  is  then  of  necessity  less  than  h  .  Let  the  lower  end  of  the  tube  be  at  a 
height  b  above  the  lower  plate. 

The  tube  is  a  sufficiently  simple  model  of  a  missile  or  aircraft  structure, 
while  the  infinite  plates  constitute  an  idealization  of  the  parallel-plate 
section  of  a  bounded-wave  simulator.  By  choosing  this  idealized  model,  it  is 
possible  to  perform  a  rigorous  and  exact  formulation  of  the  slmulator/test- 
object  interaction  problem  which  is  at  the  same  time  not  too  complicated  for  an 
analytical  solution. 


Geometry  of  the  Interaction  Between  a  Parallel-Plate  Simulator  and  a  Cylindrical  Teat  Object 


SECTION  II 


ELECTROMAGNETIC  POTENTIALS 

The  specific  goal  of  this  analysis  is  to  calculate  the  electric  current 
induced  on  the  cylinder  when  the  parallel-plate  waveguide  is  excited  by  an 
electromagnetic  wave.  This  electromagnetic  wave  can,  for  example,  be  a  simu¬ 
lated  EMP  in  the  form  of  the  transverse-electromagnetic  (TEM)  mode  of  the 
waveguide.  The  effect  of  the  simulator/test-object  interaction  is  estimated 
by  comparing  this  induced  current  with  its  corresponding  value  when  the  cylinder 
is  excited  in  free  space  in  the  absence  of  the  simulator. 

The  current  induced  on  the  cylinder  by  an  incident  electromagnetic  wave 
in  the  parallel-plate  waveguide  can  be  evaluated  by  solving  an  integro- 
differential  equation.  The  derivation  of  this  equation  requires  a  considerable 
amount  of  calculation.  One  first  introduces  the  electromagnetic  scalar  and 
vector  potentials  V  and  A  for  the  electromagnetic  fields  E80  and  Iisc 
scattered  off  the  cylinder  and  the  plates,  such  that 

Esc(r,t)  -  -  W(r,t)  -  A(r,t)  ,  Bsc(r,t)  -7*A(r,t)  (1) 

In  the  Lorentz  gauge  these  potentials  are  related  by  the  Lorentz  condition: 

7-A(r,t)  +  V(r,t)  -  0  (2) 

and  individually  satisfy  the  wave  equation: 

(v2  -  \  j-j^V(r,t)  -  0  ,  -  i-^A(r,t)  -  0  (3) 

At  z*0  and  h  ,  on  the  surfaces  of  the  two  parallel  plates,  the 
tangential  components  of  JSsc  and  the  normal  component  of  B80  must  vanish. 

By  writing  out  equations  (1)  and  (2)  in  rectangular  component  form,  one  can 
easily  convince  oneself  that  these  boundary  conditions  on  the  fields  imply 
the  following  boundary  conditions  on  the  potentials: 
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(4) 


V(r,t)  *  A  Cr.t)  "A  (r,t)  *  0 
x  ~  y 


3_ 

3z 


Az(r,t) 


0 


for  z  « 0  and  h  .  Thus  the  scalar  potential  and  the  tangential  components  of 
the  vector  potential  satisfy  the  homogeneous  Dirichlet  boundary  condition  on 
the  plates,  while  the  normal  component  of  the  vector  potential  satisfies  the 
homogeneous  Neumann  boundary  condition. 

It  will  be  convenient  to  go  over  from  the  time  variable  t  to  a  complex 
frequency  variable  s  by  introducing  the  Laplace  transforms  of  all  time- 
dependent  quantities.  For  example,  the  Laplace  transform  of  the  scalar 
potential  is 

V(_r, s)  -  j”  dt  e~StV(r,t)  (5) 

— 00 

Under  this  transformation  equations  (3)  become 

(V2- y2)V(r,s)  -  0  ,  (V2-Y2)A(r,s)  -  0  (6) 

where 

Y  -  f  (7) 

Their  solutions  can  be  represented  as  definite  integrals  over  the  surface  of 
the  cylinder: 


V(r,s)  _~l  Gjj(r,r' ,s)o(r_'  ,s)dS' 
o  * 

Ax(r,s)  -  |  GD(r,r',s)Kx(r' ,s)dS* 

Ay(r,s)  *  VQ  |  GD(r,r' ,s)Ky(r' ,s)dS' 


Az(r,s) 


GN(r,r’ ,s)Kz(r’ ,s)dS’ 


(8) 
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In  these  expressions,  a  and  K  are  the  surface  charge  and  current  densities 
induced  on  the  cylinder  by  the  incident  wave;  G^  and  are  the  Green's 

functions  for  the  Dirichlet  and  Neumann  boundary-value  problems  of  the  parallel- 
plate  waveguide;  and  the  field  position  vector  £  ■  (x,y,z)  stands  for  a  point 
in  space  while  the  source  position  vector  r_'  *  (x',y’,z')  stands  for  a  point 
on  the  surface  of  the  cylinder. 

The  Green's  functions  G^  and  G^  are  solutions  of  the  nonhomogeneous 
Helmholtz  equation: 


(V2  -  y2)GD(r,£' ,s)  -  -  53(_r-r') 

(V2-Y2)GM(r,r',s)  «  -  53(r-r') 
N  — '  — 


(9) 


They  satisfy,  respectively,  the  homogeneous  Dirichlet  boundary  condition  G^ =  0 
and  the  homogeneous  Neumann  boundary  condition  3G^/3z  *  0  on  the  plates  z=0 
and  h  .  In  addition  they  both  satisfy  the  outgoing-wave  condition  at  infinity. 
Using  standard  Green's  function  techniques,  one  finds  that 


G  (r,r’,s)  -  l  U  (x,y,x',y',s)sin  p  z  sin  p  z* 
u -  ,  n  n  n 


n“l 


with 


and 


(^(r.r'.s)  -  l  U  (x,y,x',y',s)cos  pflZ  cos  prz' 
n=0 


nir 

?n  "T 


Un(x,y,x',y',s)  *  \(  ^  +  P*  R) 


Here  K  is  a  modified  Bessel  function, 
o 


/(x-x’)2  +  (y-y')2 


and  e  is  defined  by: 
n 


2  n  ■  0 

1  n  ■  1,2,3, . . . 


(10) 

(11) 

(12) 

(13) 

(14) 
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Furthermore,  to  satisfy  the  outgoing-wave  condition,  the  branch  of  the  square 
root  in  equation  (12)  is  to  be  chosen  to  make 


as  y  -*■ «  . 


(15) 
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SECTION  III 


DERIVATION  OF  INTEGRODIFFERENTIAL  EQUATION 

The  surface  charge  density  a  can  be  eliminated  from  the  calculation  by 
using  the  continuity  equation: 

scr(_r,s)  +  V-K(r,s)  *  0  (16) 

The  independent  unknown  quantities  of  the  problem  are  then  the  two  components 
of  the  surface  current  density  R  -  One  can  set  up  two  integrodifferential 
equations  for  K  by  applying  the  boundary  conditions  on  the  tangential  compo¬ 
nents  of  the  total  electric  field  on  the  surface  of  the  cylinder.  In  a  standard 
cylindrical  coordinate  system  (p,cp,z)  with  polar  axis  coinciding  with  the  axis 
of  the  cylinder,  these  boundary  conditions  are 


E*nc(r,s)  +  E*c(r,s)  =  0  , 


Einc(r,s)  +  ESC(r,s) 
<P  ”  cp  - 


(17) 


for  _r  on  the  cylinder. 

The  two  integrodifferential  equations  so  derived  are  generally  coupled. 
However,  for  the  specific  geometry  of  this  problem,  it  is  possible  to  obtain 
from  them  a  single  closed  integrodifferential  equation  for  the  total  axial 
current  I  defined  as 


I(z,s) 


dtp  Kz(cp,z,s) 


(18) 


The  total  axial  current  by  itself  is  often  sufficient  to  characterize  the 
excitation  of  a  thin  cylinder  at  low  frequencies,  since  in  this  limit  the 
circumferential  current  is  small  by  comparison. 

To  derive  the  integrodifferential  equation  for  the  total  axial  current  I  , 
one  writes  down  the  first  of  the  boundary  conditions  (17)  in  terms  of  the 
potentials : 

V(r,s)  +  s  A  (r,s)  -  E*nc(r,s)  (19) 
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Substituting  the  integral  representations  (8)  of  the  potentials  into  equation 
(19)  and  eliminating  the  surface  charge  density  a  by  use  of  equation  (16) ,  one 
obtains 


h  f  gd(— ' * s)  [il^  Vi’’s)  +lr-Kz(i,*s>] dS< 


(20) 


-  Y 


inc . 


G  (r,r',s)K  (r',s)dS'  -  -  se  E  (r,s) 
z  ~  O  Z 


This  equation  holds  for  b  <  z,z’  <  L  +  b  ,  with 


and 


r  -  (a,<P,z)  , 


r’  -  (a,9',z') 


(•L  +  b 

dz’ 

^b 


(21) 


(22) 


It  is  to  be  especially  noted  that,  in  the  cylindrical  coordinates,  the  two 
Green’s  functions  G^  and  G^  depend  on  the  two  azimuthal  angles  9  and  9' 
only  through  the  cosine  of  their  difference  9-9*  ,  so  that 


GD(r,£',s)  -  6^(003(9-9'),  z.z’.s) 
GN(r,r',s)  =  G^(cos (9  -9' ) ,  z,z',s) 


(23) 


This  property  is  of  vital  importance  in  the  following  derivation  of  the  integro- 
differential  equation  for  I  . 

If  one  integrates  both  sides  of  equation  (20)  over  9  from  0  to  2ir  , 
two  important  events  take  place.  First,  the  9 '-dependence  of  the  Green’s 
functions  drops  out: 

•2ir 

d®  GD(cos(9 -9'),  z , z ' , s )  - 

J0 

-  2 it  rD(z,z*  ,s) 


2ir 


dif»  GD(cos  z , z ' , s) 
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j  d<P  (cos  (cp  -9*  ) ,  z,z',s)  =•  dip  GN(cos  ip ,  z,z',s) 

0  '  •'o 

-  2ir  rN(z,z',s)  (24) 


Therefore  the  subsequent  cp'-lntegration  in  equation  (20)  is  carried  out  solely 
on  the  surface  current  density  components  and  Kz  .  Second,  the  term  in 

K.  vanishes  after  the  9 '-integration  because 


9'  ■  2ir 

K  (_r '  ,s)  *  0 

?  ?'  »  0 


(25) 


On  the  other  hand  the  ^'“integration  of  Kz  yields  directly  the  total  axial 
current  I  ,  as  defined  in  equation  (18).  Gathering  together  all  these  results, 
one  obtains  the  following  nonhomogeneous  integro-differential  equation  for  I  : 


3_ 

3z 


r  L  +  b 


dz'rD(z,z' ,s)  jp-I(z',s) 


y  rL  +  b 

y  J  dz,rN(z,z',s)I(z’,s)  -  -  seqEz  (z,s) 
b 


where 


Elnc(z 

z 


1  f21r 

>3)  "  27 


d<P  Eznc(a,9,z,s) 


(26) 


(27) 


That  is,  E  is  the  average  value  of  the  z-component  of  the  incident  electric 
z 

field  around  the  circumference  of  the  cylinder  at  a  fixed  value  of  z  . 

One  should  emphasize  that  equation  (26)  is  an  exact  equation.  No  approxi¬ 
mation  has  been  invoked  in  its  derivation.  This  remarkable,  rigorous  result 
is  actually  made  possible  through  the  judicious  choice  of  the  geometry  of  the 
model.  If  a  different  geometry  vere  chosen,  it  is  doubtful  that  a  closed 
equation  for  I  similar  to  equation  (26)  could  be  derived  without  the  help 
of  assumptions  and  approximations. 
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The  two  functions  and  rN  in  equation  (26)  can  be  evaluated 
explicitly  in  closed  form.  Substituting  equation  (10)  into  equation  (24) , 
one  obtains 


Use  has  been  made  of  the  addition  theorem  for  the  cylindrical  functions  (ref.  2) 
in  evaluating  the  integral  in  equation  (29) . 


SECTION  IV 


REDUCTION  TO  MATRIX  EQUATION 

Equation  (26)  is  a  second-order  integrodif ferential  equation.  It 
therefore  requires  the  specification  of  two  boundary  conditions  in  the 
coordinate  z  for  the  unique  determination  of  the  total  axial  current  I  . 
These  boundary  conditions  can  be  taken  as  the  two  end  conditions. 

I(z,s)  -  0  for  z-b  and  L  +  b  (30) 

However,  it  is  often  inconvenient  to  treat  the  problem  in  hand  directly  as  a 
two-point  boundary-value  problem  for  equation  (26) .  For  the  purposes  of 
analytical  approximation  or  numerical  solution,  much  advantage  can  be  gained 
by  combining  equation  (26)  and  the  boundary  conditions  (30)  and  reducing  them 
to  a  single  matrix  equation. 

The  transition  to  a  matrix  equation  is  achieved  by  going  over  to  a 
discrete  system  of  description.  To  this  end  one  Introduces  a  complete  set 
of  mutually  orthogonal  functions  over  the  Interval  b  <  z  <  L  +  b  which 
vanish  at  the  two  end  points  z-b  and  L  +  b  .  This  set  consists  of  the 

sine  functions  sin  k  (z-b)  with 

m 


k 


m 


m  -  1,2,3, .. . 


(31) 


These  functions  can  be  considered  as  the  orthogonal  basis  vectors  of  an  infinite¬ 
dimensional  vector  space.  The  current  I  ,  satisfying  the  boundary  conditions 
(30),  can  be  expanded  in  terms  of  these  basis  vectors: 

OB 

I(z,s)  -  7  C  (s)sln  k  (z-b)  (32) 

,  n  ® 

m-l 

The  discrete  set  of  expansion  coefficients  C  completely  determines  I  . 

m 

Next,  one  substitutes  equations  (28)  into  equation  (26)  and  obtains 
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(33) 


where 


I7  l  W  (s)A  (s)sin  p  z  ~  V  l  W  (g)B  (s)cos  p  z 
3*  nt1  n  n  n  ^0  n  n  n 


n  n-0  n  n 


A  (s)  •  dz'sln  p  z'  |-t  I(z',s) 

n  dz 


rL 

i  (•>  • 
n 


dz’coe  p^z'Kz'  ,s) 


Performing  an  integration  by  parta  in  equation  (34)  and  making  use  of  equation 

(30),  one  finds  that  A  and  B  are  related: 

a  n 

An(8)  *  “  PnBn(9)  (36) 

Substituting  equation  (32)  into  equation  (35),  one  has 


’n(,)  '  J  “a,'.'*’ 


where 


rL  +  b 

dz'cos  p  z'  sin  k  (z'  -b) 
)  n  m 


Note  that  the  index  n  runs  from  0  to  »  while  the  index  m  runs  from 
1  to  «  . 

Equation  (33)  can  be  regarded  as  a  vector  relation  in  the  infinite- 
dimensional  vector  space.  The  component  of  this  relation  along  the  basis 

vector  sin  k  (z  -  b)  can  be  obtained  by  multiplying  equation  (33)  by 

m 

sin  k  (z-b)  and  integrating  over  z  from  b  to  L  +  b  .  The  derivative 

n 

in  the  first  term  is  again  treated  by  integration  by  parts.  The  result  then  is 


l  Bnmk«Wn(8)An(s>  +  ^  l  °nmV8)Bn(8)  “  Dm(s) 
n“l  n*0 


Bin' - 
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SECTION  V 


APPROXIMATE  ANALYTICAL  SOLUTION 


The  matrix  equation  (43)  is  a  concise  and  exact  formulation  of  the 
simulator/test-object  interaction  problem.  From  the  standpoint  of  solving 
the  problem  it  is  also  a  convenient  starting  point.  For  example,  an  accurate 
numerical  solution  can  always  be  obtained  by  using  standard  methods  in  numerical 
matrix  calculus.  In  the  following,  however,  an  approximate  analytical  solution 
will  be  constructed  instead. 

The  approximation  employed  below  is  based  on  the  observation  that  the 
diagonal  elements  of  the  matrix  M  in  equation  (43)  are  considerably  greater 
than  the  off-diagonal  elements.  The  reason  behind  this  disparity  is  as  follows. 
If  one  examines  closely  the  structure  of  the  matrix  element  M  given  in 

01  x, 

equation  (44) ,  one  finds  that  the  combination  a  a  has  the  most  violent 

mo  nx 

variation  with  the  summation  index  n  .  It  has  therefore  a  great  effect  on 
the  value  of  M  .  The  factor  a  can  be  evaluated  from  equation  (38): 


a 

run 


£cos  pnb  - 


(-l)mcos  pn(L  +  b)j 


(45) 


It  assumes  both  positive  and  negative  values.  In  an  off-diagonal  matrix  element 

M  .  with  m  ¥  i  ,  these  positive  and  negative  values  of  a  and  a  .  to  a 
cl  r  run  nH 

considerable  extent  mutually  cancel  out  during  the  summation  over  n  .  However, 
the  situation  is  different  with  a  diagonal  matrix  element  M  .  With  m  *  l  , 

2 

the  product  a  a  .  becomes  the  square  a  which  is  always  positive.  It 
nm  nx  nn 

contributes  constructively  to  M  ,  and  adds  up  to  a  large  value.  For  small 
values  of  m  and  1  ,  sample  numerical  calculations  show  that  the  diagonal 
elements  are  one  order  of  magnitude  greater  than  the  off-diagonal  ones.  The 
disparity  is  expected  to  increase  for  larger  values  of  m  and  l  . 

This  property  of  the  matrix  elements  is  the  basis  of  an  approximation 
scheme  for  solving  the  matrix  equation  (43) .  One  first  rearranges  the  equation 
in  the  fora 

00 

Mma(s)Cm(s)  -  DB(s)  -  l  MnJ(s)C]l(s)  (46) 

Jt-1 

fc#m 
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The  diagonal  element  has  been  taken  out  from  the  sum.  The  remaining 
sum  contains  only  off-diagonal  elements,  and  is  regarded  as  a  perturbation. 
A  first  approximation  to  the  solution  of  equation  (46)  is  therefore 


M 


D  (s) 

IB 

TT> 


(47) 


Successive  approximations  can  be  obtained  by  Iteration.  For  example,  the 
second  approximation  is 


M 


D  (a) 
tn 

TZY 


mm 


I 

4- 1 
i^m 


V8> 

H~_(b) 


(48) 


In  the  following  sections  the  solution  in  the  first  approximation  will  be 
examined . 


SECTION  VI 


COMPLEX-FREQUENCY-PLANE  SINGULARITIES 


In  the  first  approximation  of  equation  (47)  the  total  axial  current 
induced  on  the  cylinder  is  given  by 

I(z,s>  ’  l  STTsT  sin  km(z  ~ b)  (49) 

m«l  mm 

In  the  time  domain  this  current  becomes 


I(z,t) 


•C+j« 

C-j« 


of  £  0  (s) 

dS  e  i  H~U)a±n  km(z-b) 

m-1  mm 


(50) 


where  C  is  an  appropriate  positive  constant.  The  value  of  the  integral 
is  dependent  on  the  singularities  of  the  integrand  in  the  left  half  of  the 
complex  s-plane.  The  parameters  of  these  singularities  can  be  used  to 
characterize  the  excitation  of  the  cylinder. 

The  s-plane  singularities  can  be  grouped  into  three  classes.  The  first 

class  consists  of  singularities  of  the  factor  D  .  By  equation  (41)  these 

m 

singularities  are  introduced  by  the  frequency  dependence  of  the  Incident 
luc 

electric  field  Ez  .  The  second  class  consists  of  singularities  of  the 
factor  M  .  By  equation  (44)  these  singularities  are  contained  in  the  func- 
tion  defined  in  equation  (29),  From  the  well-known  analytical  properties 

of  the  modified  Bessel  functions,  one  concludes  that  these  singularities  are 
branch  cuts  with  branch  points  determined  by  the  vanishing  of  the  argument  of 
Che  function  KQ(/f^  +  p^  a)  .  By  equation  (7)  the  branch  points  are  located  at 

s  -  ±  J  Sjp  n  -  0,1,2,...  (51) 

along  the  imaginary  axis.  They  are  characteristic  of  the  parallel-plate 
waveguide,  and  depend  neither  on  the  length  L  nor  the  radius  a  of  the 
cylinder.  In  fact  they  correspond  to  the  excitation  thresholds  of  the  wave¬ 
guide  propagation  modes.  Examples  of  these  waveguide  mode  singularities  are 
shown  in  figure  2. 
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Figure  2.  Branch-Cut  Singularities  in  the  Complex-Frequency  Plane 
Characterizing  the  Waveguide  Modes 


1  ^ 


The  third  class  of  s-plane  singularities  consists  of  the  zeros  of  M 

mm 

They  give  rise  to  poles  in  the  Integrand  of  equation  (50),  and  are  determined 
by  solving  the  transcendental  equation 

M™<9>  •  I„(4+  ‘  0  ■■  1.2.3. •••  <52) 

n»0  \  c  / 


The  solutions  depend  on  all  the  parameters  h,  L  and  a  of  the  model.  They 
can  be  interpreted  as  the  natural  frequencies  of  certain  natural  modes  of  the 
surface  current  on  the  cylinder.  For  each  m  ,  there  are  two  solutions  of 
equation  (52)  which  form  a  complex-conjugate  pair.  They  lie  in  the  second 
and  third  quadrants  of  the  s-plane.  The  natural  modes  corresponding  to  these 
natural  frequencies  are  the  sine  functions  in  equation  (50).  Thus,  in  this 
approximation,  the  natural  modes  are  simple  sinusoidal  currents. 

Let  the  solution  of  equation  (52)  in  the  second  quadrant  be  denoted  by 

* 

s  .  The  other  solution  is  its  complex  conjugate  s  in  the  third  quadrant, 
m  m 

Decompose  into  its  real  and  Imaginary  parts: 


m 


m 


+  i  8m 
m 


(53) 


In  principle,  s^  and  s”  can  be  obtained  accurately  by  numerical  techniques. 
In  the  following,  however,  one  will  be  content  with  deriving  approximate 
analytical  expressions.  The  method  of  approximation  will  be  based  on  the 
assumption  that  the  two  poles  lie  close  to  the  imaginary  axis,  so  that  they 
can  be  approximately  located  by  observing  their  influence  on  the  behavior  of 
Mmm  along  the  imaginary  axis. 

On  the  imaginary  s-axis,  let 

s  -  ju  ,  y  ■  jk  (54) 

where  u  and  k  are  both  real  and  vary  from  -“  to  «  .  They  are  related  by 

u  ■  kc  (55) 
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The  function  W  can  be  decomposed  into  its  real  and  imaginary  parts: 
n  - 


W  (jco)  -  W'(u>)  +  j  Wm(uj) 
n  n  n 


By  equations  (15)  and  (29) ,  one  has 


w'(co)  -  rrr-  i 


(to)  *  — r —  I  \  *p  -  k  a/iC\rp  -  k  a  /  0  l  p  -  k 
n  irhe  ov  rn  'o'  *11  '  v  Kn 

n 


)  Ki^n-k2  a)e(p^-k2 


-i-  J  (»(2-  p2  a 
2ne„  o  *n 


*2  r  ),f,2  2) 

c  -Pn  a /e\k  -pn/ 


W"(ai) 

n 


w  1  t2//2  2  \ft(,2  2\ 
; — -7 —  J\rk  -p  a  /  0  \  k  -p  / 
M2hen  °'  n  '  n  ’ 


where  J  and  Y  are  the  Bessel  and  Neumann  functions,  and  0  is  the  unit 
o  o  ’ 

step  function: 

(l  x  >  0 


0  x  <  0 


The  square  roots  in  equation  (57)  are  all  considered  positive. 

In  a  similar  manner  the  value  of  M  along  the  imaginary  axis  can  be 

mm 

decomposed  into  its  real  and  imaginary  parts: 


M  (ju)  =  M»  (co)  +  j  M"  (to) 
nun  Hull  nun 


M'  (w)  *  \  ( -  ^7  + 

mm  n”0  \  c  n/  ”  M 


M"  (w)  =  I  I  -  +  Pn  )W"(w)a 

111111  n-0  \  CT  n/  "  " 


When  j«  is  close  to  the  solution  s^  of  equation  (52)  in  the  second  quadrant, 
can  be  seen  to  pass  through  a  zero.  This  property  is  used  to  determine 
s”  approximately: 
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which  is  a  transcendental  equation  involving  real  quantities  only.  It  is 
simpler  to  solve  than  the  original  complex  equation  (52). 


A  practical  technique  to  solve  equation  (61)  is  by  iteration.  Using 
equation  (60)  one  can  rewrite  it  as 


_»f  _  _ 

s  -  c 
m 


n*0 


2  2 
p  W'(s")a 
n  n  m  nm 


W1  (sf,)a 
n  tn  nm 


'n*0 


(62) 


The  iteration  is  performed  by  successively  substituting  trial  solutions  into 
the  right-hand  side  and  deriving  improved  approximations  on  the  left-hand  side. 
In  the  appendix  it  is  shown  that,  in  the  infinitely  thin  cylinder  limit,  the 
solution  of  equation  (62)  is  simply 


for  f  ■* 

L 


(63) 


Mence,  for  a  thin  cylinder,  one  can  use  this  result  to  start  the  iteration. 
An  improved  approximation  to  s"  is  therefore 


s 


M 

m 


(64) 


The  following  calculation  will  be  based  on  this  approximate  solution  of 
equation  (61) . 

In  the  neighborhood  of  this  solution,  M'  has  the  following  approxi- 

mm 

mate  expansion: 


where  is  given  by  equation  (64)  and 


is  a  real  quantity  defined  by 


Q 


m 


00 


l 

n"0 


W'(k  c)a 
n  m 


2 

ntn 


(66) 


The  imaginary  part  has  a  hump  around  this  solution  and  can  be  approximated 

by  its  value  at  u>  =  k  c  : 

m 


M"  (m) 
mm 


M"  (k  c) 
mm  m 


(67) 


Gathering  together  these  results,  one  finds  that  the  complex  function  M  has 

TTWT1 

the  expansion 


where 


Mmm(ju)>  = 
mm 


ft  (D)-s"  +  js') 
m  m  J  m 


for 


m 


M"  (k  c) 
m  s«  mm  m 
m 


(68) 


(69) 


Equation  (68)  shows  that,  in  the  complex  s-plane  (jm  •*  s)  ,  the  function 

M  has  a  zero  at  s  ■  s'  +  1  s"  ,  Therefore  the  reciprocal  of  M  has  a 

mm  m  m  m  mm 

pole  at  the  same  point.  In  the  neighborhood  of  this  point,  one  has 


-j 


M  (s)  "  n  (s-s  ) 
mm  m  m 


for  s  =  s 


m 


(70) 


Similarly,  by  repeating  the  calculation  in  the  third  quadrant,  one  finds  that 

M  has  a  second  zero  at  s  ,  and  that 
mm  m  ’ 


1 


M 

mm 


(s) 


Qjs  ~  O 

m  m 


for  s-s* 
m 


(71) 


In  summary,  the  total  axial  current  on  the  cylinder  given  in  equation 
(49)  has  the  following  expansion  in  the  complex  s-plane: 
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+  f(z,s) 


(72) 


I(z,s) 


l 

m*l 


R  (z)  R*(z) 

m  in 

s  -  s  ~ 


m 


s  -  8 


m  J 


It  consists  of  a  sum  of  simple  poles  and  a  remainder  term  V  .  For  each 
value  of  m  there  is  a  pair  of  poles  situated  at 


s  *  s'  +  ]  s" 
in  m  m 


and 


* 

s 

m 


s’ 

m 


3  3" 

in 


(73) 


where  and  s”  are  approximately  given  by  equations  (64)  and  (69).  By 

eauations  (41),  (49)  and  (70),  the  residue  R  is 

tn 


Rmu) 

m 


3  D  (s_) 


0 


m 


sin  k  (z  -  b) 
m 


Js  e 
m  o 


m 


sin  k  (z  -  b) 
m 


r 


dz’sin  k  (z*  -  b)Einc(z' ,s  )  (74) 

m  z  m 


These  poles  correspond  to  the  natural  modes  of  the  cylinder.  The  remainder 
V  contains  all  other  contributions  from  the  branch  cuts  of  the  waveguide 
nodes  as  well  as  whatever  singularities  the  incident  electric  field  may 
introduce. 
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SECTION  VII 


NUMERICAL  RESULTS 


The  three  quantities  s'  ,  s'*  and  ft  ,  characteristic  of  the  simple-pole 

m  m  m 

singularities  in  the  complex  s-plane,  are  given  analytically  in  equations  (64) , 
(66)  and  (69).  In  the  following  they  are  evaluated  numerically  as  functions 
of  the  geometrical  parameters  h  ,  a  ,  L  and  the  natural-mode  index  m  .  One 
takes  a  case  in  which  the  cylinder  is  situated  midway  between  the  two  parallel 
plates,  so  that 

b  *  (h-L)/2  and  L  +  b  -  (h  +  L)/2  (75) 


This  is  a  highly  symmetric  situation.  Accordingly,  the  quantity  a  assumes 

nm 

a  simple  form: 


2k 


m 


nm  2  2 

k  -  p 
m  rn 


p  L 

r  ntr  *n 

'cos  -j-  cos  -y-  m  odd,  n  even 


P  L 

.  nir  .  rn 

sin  *Tjj-  sin  m  even,  n  odd 


(76) 


*2 


otherwise 


Thus  is  nonzero  only  if  n  and  m  are  of  opposite  parities. 

The  numerical  calculation  is  performed  for  the  first  six  natural  modes 
(m  ■  1  to  6)  .  One  picks  four  values  of  the  ratio  a/L  : 


jr  -  0.003,  0.01,  0.03,  0.1  (77) 

and  three  values  of  the  ratio  L/h  : 

£-  0.2,  0.5003,  0.8003  (78) 

The  results  are  shown  in  tables  1  to  3.  They  are  also  plotted  in  figures  3  to  5 

The  odd  choices  of  the  last  two  values  of  L/h  in  equation  (78)  need 
immediate  explanation.  When  the  ratio  L/h  is  given  simple,  rational  values 
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Table  1 


POLE  SINGULARITY  PARAMETERS  FOR  L/h-0.2 


26 


Table  3 


*  ’4fc  4’ 


such  as  0.5,  it  may  happen  that,  for  a  certain  pair  of  integers  m  and  n 
the  following  equality  holds: 


K  '  pn 


(79) 


In  this  event  one  finds  a  logarithmic  divergence  in  equation  (64)  arising  from 

the  exact  vanishing  of  the  argument  of  the  function  Kq  or  Yq  .  This 

divergence  nevertheless  is  not  really  there,  and  can  be  traced  back  to  a 

harmless  term  vanishing  as  e  2-n  e  in  the  function  M  .  It  is  introduced 

mm 

into  the  picture  through  the  iterative  scheme  of  equation  (62)  and  the  use  of 

k  c  as  the  first  trial  solution  in  equation  (64) .  It  can  be  easily  avoided 
m 

by  starting  the  iteration  with  a  slightly  different  trial  solution.  An  alter¬ 
native  procedure  is  to  alter  the  ratio  L/h  by  a  small  amount.  For  example, 
the  ratio  L/h  *0.5  can  be  replaced  by  L/h  *  0.5003  ,  thereby  destroying 
the  equality  (79) . 

For  m  -  1  to  6  and  n  of  opposite  parity  to  m  ,  the  equality  is  not 
satisfied  for  L/h  *  0.2.  However,  for  L/h  ■  0.5  ,  equation  (79)  holds  for 
m  -  1,  3  and  5  .  For  L/h  *  0.8  ,  it  holds  for  m  »  4  .  Accordingly,  in  the 
numerical  calculation,  the  last  two  simple  values  of  L/h  are  changed  slightly 
to  0.5003  and  0.8003. 

Even  though  no  divergence  really  occurs,  equation  (79)  still  clearly 
represents  a  condition  of  resonance  and  should  lead  to  observable  effects. 

The  resonance  results  from  the  equality  of  the  wavelengths  of  certain  charac¬ 
teristic  excitations  on  the  cylinder  and  in  the  waveguide.  Its  effects  are 
auite  evident  in  figures  3  to  5.  At  precisely  those  values  of  m  enumerated 
in  the  preceding  paragraph,  the  complex  poles  show  a  decided  shift  toward 
the  imaginary  axis.  This  means  that  the  corresponding  natural  modes  of  the 
cylinder  are  only  weakly  damped.  Thus,  under  the  excitation  of  a  broadband 
incident  wave,  the  response  of  the  cylinder  will  be  predominantly  in  these 
special  modes. 

All  the  results  in  tables  1  to  3  are  of  course  obtained  by  using  the 
first-iterated  solution  of  equation  (62).  The  quality  of  this  solution  can 


be  measured  by  comparing  the  output  value  of  a”  with  its  input  trial  value 

k  c  .  The  tables  show  that  these  two  values  are  very  close.  Specifically, 
m 

the  calculated  value  of  s"l/itc  is  very  close  to  m  .  This  agreement  can  be 

tn 

taken  as  an  indication  of  the  reliability  of  the  first-iterated  solution. 
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SECTION  VIII 


CONCLUSIONS 

It  is  possible  to  draw  a  number  of  general  conclusions  from  the  specific 
numerical  results. 

The  effects  of  the  simulator/test-object  interaction  are  clearly  seen 
in  figures  3  to  5.  The  locations  of  the  poles  in  the  complex-frequency  plane, 
corresponding  to  the  natural  modes  of  the  current  on  the  cylinder,  are  strongly 
dependent  on  the  length-to-separation  ratio  L/h  .  When  L/h  is  small,  the 
dimensions  of  the  simulator  are  much  greater  than  those  of  the  test  object. 

The  simulator/test-object  interaction  is  weak,  at  least  for  the  lowest  few 
natural  modes.  For  example,  the  distribution  of  the  poles  in  figure  3  for 
L/h  *  0.2  is  not  substantially  different  from  that  in  the  case  of  a  thin 
cylinder  in  free  space  (ref.  3).  As  the  ratio  L/h  is  increased,  the  effects 
of  the  interaction  strengthen:  the  poles  are  shifted  violently  around.  The 
dependence  of  their  movement  on  L/h  ,  however,  is  not  monotonic  but  rather 
oscillatory  and  in  the  nature  of  resonances.  For  a  given  ratio  L/h,  modes 
that  satisfy  the  resonance  condition  (79)  are  heavily  favored  by  the  simulator/ 
test-object  interaction,  in  the  sense  that  they  are  only  weakly  damped.  These 
resonant  cases  are  exemplified  by  the  numerical  results  in  figures  4  and  5. 

The  poles,  however,  are  not  sufficient  in  themselves  to  provide  a  complete 
measure  of  the  simulator/ test-object  interaction.  The  presence  of  the  simulator 
structure  indroduces  additional  singularities  into  the  complex- frequency  plane. 

For  the  parallel-plate  bounded-wave  simulator  in  this  problem,  these  singularities 
consist  of  an  infinite  number  of  branch  cuts.  The  contribution  from  the  branch 
cuts  must  be  considered  if  one  is  to  obtain  a  full  characterization  of  the  simu¬ 
lator/test-object  interaction. 

The  numerical  results  of  this  study  are  based  on  an  approximate,  analytical 
formula  for  the  solution  of  an  exact  matrix  equation.  It  will  be  worthwhile 
to  attempt  an  accurate,  numerical  solution  of  the  matrix  equation  directly. 

Such  a  solution  is  feasible  with  the  availability  of  large,  advanced  computers. 

It  will  not  only  provide  accurate  technical  data,  but  will  also  represent  a 
valuable  standard  whereby  the  reliability  of  future  analytical  investigations 
can  be  gauged. 
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APPENDIX 


THIN-CYLINDER  LIMIT 

In  the  thin-cylinder  limit  (a/L  •+  0)  equation  (62)  can  be  shown  to  have 

the  simple  solution  s"  =  i  k  c  . 

mm 

2  -4 

By  equation  (45)  the  quantity  varies  like  n  for  large  n  .  Thus 

the  contribution  to  the  two  sums  in  equation  (62)  comes  mainly  from  terms  of 
small  n  .  Because  of  this  property,  one  can  replace  the  function  by  its 

thin-cylinder  limit  at  small  n  .  By  equation  (57)  one 

sr- *»(!)  as  !* 

n 

Therefore  equation  (62)  becomes 
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The  two  sums  can  be  evaluated  by  invoking  equations  (38)  and  (40) : 
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By  equation  (42)  it  can  also  be  rewritten  as 
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where  use  has  been  made  of  the  identities 
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for  b  <  z,z'  <  L  +  b  .  Substituting  equations  (A4)  into  equation 
obtains  the  simple  result 
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